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Abstract 

We obtain the explicit expression of the friction on a moving quan- 
tized vortex, following the formulation of Thouless, Ao and Niu. It is 
shown that the friction on a moving vortex is sensitive to details but 
does not change the transverse force. We provide a general thermo- 
dynamic interpretation for the mutual independence of the transverse 
force and the friction. The friction is evaluated for the case of quasi- 
particle contributions in a clean fermionic superfluid, showing a new 
feature of logarithmic divergence. 
PACS#s: 67.40.Vs; 47.37.+8 

In physical systems few qualities have generic values independent of ma- 
terial details. When such cases occur, usually there are either symmetry or 
topology reasons]]]]. The controversy over dynamics of quantized vortices 
in superfluids or superconductors 0, || may be partly understood as due to 
the apparent conflict between the exact results @, || § with the existence of 
various scattering processes in the system. When a vortex moves, in addition 
to the hydrodynamic Magnus force which is proportional to the superfluid 
density, the coupling between the vortex to quasiparticles, phonons and im- 
purities may cause extra forces in both transverse and longitudinal directions 
of the vortex motion. Those possible non-vanishing forces have always been 
the underlying assumption for phenomenological models used in analyzing 
experiments]^]. |7|. 
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While the transverse force on a moving vortex is controversial, the fric- 
tion is poorly studied. The usual quoted friction formula derived from a 
microscopic Hamiltonian is obtained through a relaxation time approxima- 
tion in force-force correlation functions. However, this type of procedure has 
been known to be wrong since 60's, in the context of obtaining the friction 
of electrons, i.e. resistivity. The lack of explicit and correct microscopic 
derivation of friction has two consequnces. First, though the friction of vor- 
tex motion is an important quantity in many experiments, it has not been 
received its due attention. Second, the incomplete understanding of the fric- 
tion has generated doubts on the exact result on the transverse force. 

A self-contained theory of both the friction and the transverse force by a 
rigorous and elementary method is undoubtedly desirable. To provide such 
a calculation in an explicit and complete manner is the main purpose of the 
present letter. The method we choose follows that of Ref. B. In order 
to obtain the friction, we use a limiting process similar to the one in the 
transport theory and in non-equilibiurm statistical mechanics. This limiting 
process has only been mentioned but not explicitly discussed in Ref. due 
to its sensitivity to details. 

In the following, after giving a general expression for the friction, as an 
example, we evaluate the friction for the case of a clean fermionic superfluid. 
The resulting friction is completely new. It comes from the off-diagonal 
potential scattering of the extended quasiparticles, and is stronger than the 
ohmic damping by a logarithmic diverging factor. This new result illustrates 
the directness and usefulness of our formulation, another purpose of the 
present paper. 

We consider an isolated vortex in the superfluid whose position ro is 
specified by a pinning potential. The system is otherwise homogeneous and 
infinite. There is no externally applied supercurrent or normal current. The 
vortex is allowed to move slowly. Hence the system Hamiltonian H con- 
tains a slow varying parameter r (t). The many body wavefunction of the 
superfluid \^f a (t)) can be expanded in terms of the instantaneous eigenval- 
ues E a (r Q ) and eigenstates |^ a (ro)), for which we choose phases such that 
(4>a\ipa) = 0. Because of our assumption of a homogeneous superfluid, E a (r ) 
is independent of both r and time t. With those considerations, the many- 
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body wavefunction \ty a (t)) can be expressed by 

|*«(f)> = e-^^|^(ro)) + E MOe^^llMro)) , (1) 

where, to the first order in velocity, a a (t) = 1, and 

ac{t) = - fdt\M^)j {Ea, - Ea)t ' /h ■ (2) 
Jo 

Here V denotes the partial derivative with respect to the position r of the 
pinning potential, the vortex position. This gives the expectation value of 
the force on the vortex as 

F = -£/ Q <tt Q |V tf|tt«) 

a 

= -£/«^«lVoif|^> 



L 



E /a(^(roW)|V //|Va'(ro(i))) x 

*^'(^(r (0)lV'a(ro(i')))e i(E «'- E " )( *'-* )/7i + c.c., (3) 



where / a is the occupation probability of the state a. For a vortex moving 
with a small and uniform velocity Vy, 

|^a(r )) = v v • |V ^a(r )) , (4) 

The first term in the right hand side of Eq.(3) is independent of vy and will 
be ignored. In fact, it is zero because of the translational invariance regarding 
to the vortex position. The integration in time can be carried out directly. 
In addition, we also need to use the following relations, 

(V>a|V #|Vv> = (E a , - E a ){i) a \V^ a >) = (E a - E a ,){VMi) a .) , (5) 

which are obtained by taking gradient Vo with respect to H\ip a i) = E a i\ip a i) 
and (ip a \H = E a (ip a \, then multiplying from left /right by (ip a \ or \ip a ') re- 
spectively After making simplifications by using Eq.(5), Eq.(3) gives 

^ 1 v-^ fa -fa' sin(E a - E a >)t/h 2 

F = 9 £ p p (p p Wfi \M^ H\M\ Vy 
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+ih ^2 fa (1 - COs{E a - E a ,)t/h) {((^ q |V t/v) X (Vo^'a'lV'a)) ' z} Vy X Z , 

(6) 

the longitudinal and transverse forces on a moving vortex. 

To obtain the long time behavior of Eq.(6) requires a limiting procedure. 
There are two ways of taking limiting sequences. A possible sequence is 
to take the low frequency limit before the thermodynamic limit. As we 
will find out, the transverse force is independent of the limiting process. 
Therefore such a limiting sequence is correct as long as only the transverse 
force is concerned. In such a calculation, because the energy levels have been 
treated as discrete ones, there is no friction. In order to obtain friction, we 
have to take the thermodynamic limit before the low frequency limit, the 
correct limiting procedure. For this purpose, we use the Laplace average 0, 
lim^oo F(t) = lim e ^ 0+ F(e), where F(e) = e/ °° dt F(t)e~ et . We will also use 
the identity l/(e + ix) = tc5(x) — iP(l/x). For the transverse force Fj_, we 
have 

poo 

lim e / dte- et (1 - cos(£ Q - E a ,)t/h) = lim (1 - n5(E a - E a ,)e) = 1 , 

e^0+ JO e^0+ 

because the summation over states in Eq. (6) is well behaved regardless whether 
E a is discrete or continuous. Hence 

F±=iH^2 f a {((V^|VoVv) X (V 1pa'\lpa)) ' z} Vy X Z . (7) 

a'j^a 

This is precisely what has been obtained in Ref . , where further calculations 
lead to F± = —hLp s vv x z, independent of details of the system. Here p s is 
the superfluid number density, L the length of the vortex. The longitudinal 
force, friction, is given by Fii = — rjvy with 

V = ^J2 n l a ~J^ 5(E a -E af )\(rp a \V H\M\ 2 . (8) 

^ a'^a a ' a 

The coefficient of friction rj is determined by low energy excitations such 
as phonons, extended quasiparticles, and localized quasiparticles when their 
discrete energy spectrum is smeared out by impurities. This expression is 
identical to the result in Ref. |1 for the case of ohmic damping in the zero 
frequency limit. Eq.(8) will not pick up any superohmic contributions, and 
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will give infinity for any subohmic contributions. We point out that while 
with the aid of Eq.(5) the transverse force can be expressed only in terms 
of the wavefunction or the density matrix without explicitly refering to the 
Hamiltonian or its eigenvalues, as showing by Eq.(7), the longitudinal force, 
the friction, cannot, as showing by Eq.(8). The explicit dependence on the 
Hamiltonian or its eigenvalues in Eq.(8) is the source of the sensitivity of 
friction to details of the system. 

For the superfluid 4 He, there is no full microscopic theory yet [7| . We will 
not attempt to evaluation of Eq.(8) for this superfluid here because of the 
sensitivity of the friction to details. The situation can be very different in 
the case of the recent Bose-Einstein condensed systems, where a well defined 
microscopic theory is supposed to be known. Instead, as an example to 
illustrate the directness and usefulness, we evaluate Eq.(8) for the case of a 
homogeneous fermionic superfluid using BCS theory with s-wave pairing. At 
finite temperatures the extended states above Fermi level, the quasiparticles, 
are partially occupied. The vortex motion causes transitions between these 
states and gives rise to friction. The transitions between different single 
quasiparticle levels (ip a \ VoH\ip a /) are considered here since they dominate 
the low energy process. The quasiparticles are described by the eigenstates, 
u a and v a , of the Bogoliubov-de Gennes equation. Their behaviors in the 
presence of a vortex has been well studied in Ref . fT0|| . We may take 

M = ( ij^ ) = -Le^j^Wfir) , (9) 

with r measured from the vortex position, and 9 the azimuthal angle around 
the vortex. In order to obtain an analytical form for the transition element, 
we use a WKB solution for /(r), 




'e 2 - 


|A(r) 


I 2 , 


'e 2 - 


\A(r) 


I 2 , 



f(r) = x 



exp {i I dr' [ k 2 ;—^- ± |? V^ 2 - l A ( r ')l 2 ) \ + c - c - ( 10 ) 



1/2' 



Here k 2 = k 2 — k 2 , R is the radial size of the system. This WKB solution 
is valid when r is outside the classical turning point r t = \y\/k p . Here r t is 
the impact parameter. A WKB solution also exists inside the turning point. 
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However, because it approaches zero as (rk p y^/\fj,\\, the contribution to the 
transition elements from this region is small, and will be set to zero. The 
transition elements are then given by 



|(^ a |V ##o 



r 2 
dx(u* a ,(x)V Q Av a (x) + v* a ,(x)V Q A*u a (x)) 

A 4 

, > i kp 



Here A^ is the value of |A(r)| far away from the vortex core. Physically 
it means that if the classical quasiparticle trajectory is far away from the 
vortex core, it will not contribute to the friction. The summation over states 
in Eq.(8) is replaced by 



E = E / dE i dE 2 
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Ex E 2 ( m \ R 
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after considering the density of states. 

Substituting Eq.(ll) into Eq.(8), using the quasiparticle distribution func- 
tion f a = \j{eP Ea + 1), the coefficient of friction is given by 

_ Lm%At(3 f°° E 2 1 ( . 

V 8n 2 h 3 Ja E 2 - Al E 2 cosh 2 (J3E/2) ' 1 ' 

The integral in Eq.(12) diverges logarithmically. It implies that the spectral 
function corresponding to the vortex-quasiparticle coupling is not strictly 
ohmic but has an extra frequency factor proportional to ln(Aoo/^a;)|]ll], |]. 
When hu is not very small comparing to A^, which may be realized when 
close to T c , we can ignore the logarithmic divergence in Eq.(12) by using the 
density of states for normal electrons to obtain a finite friction, i.e. replacing 
E 2 j{E 2 — A 2 ^) with 1 in Eq.(12). Close to T c , the friction approaches to 
zero the same way as A 2 ^, which is proportional to the superfluid density p s . 
When — ln(fia;/A 00 ) is large, we need to use a more accurate expression of 
vortex friction obtained in Ref. ||. Straightforward evaluation shows that 
in such a case 

Lm 2 i A^l3 1 
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Here u c is the low frequency cut-off. It is determined by the size of the system 
for a single vortex, and by the inter-vortex distance for a vortex array. 

It should be emphasized that the logarithmic divergence comes from the 
interplaying between the divergence in the density of states and the off- 
diagonal potential scattering. We can consider a situation that we physically 
create a pinning center to trap the vortex and guide its motion. In such 
a case the vortex has a diagonal potential. If the scattering is dominated 
by the diagonal potential, e.g., by the pinning potential, additional factor 

coming from \u a \ 2 \ v a\ 2 will remove this logarithmic divergent. This 

again shows the sensitivity of the friction to system details. 

This extended state friction contribution exists for both clean and dirty 
superconductors at finite temperatures. Close to the transition temperature, 
it scales linearly with the superfluid density, and is exponentially small when 
T « Aqo. For intermediate temperature T ~ Aoo, using £ ~ h 2 k F /mA^ 
and iV(0) = rnkp/^fi 2 , i] ~ LhN "(0) / 'k bT '. For dirty superconductors 
additional contribution to the friction arise from vortex core states. This 
contribution, which remains finite at zero temperature, was first discussed 



phenomenologically in Ref . ||1 2|| . Eq.(12) also permits us to do a detailed 
microscopic calculation for the core contributions, as pointed out in Ref. 0. 
Another application of the above results may be the dynamics of the axisym- 
metric vortex in 3 He B phase, where the Bogoliubov-de Gennes equation is 
essentially the same as the one in a s-wave superconductor [p~3 . 



As we have shown, the friction on moving vortex has certain interesting 
features which can be tested experimentally. It is also clear that transverse 
force is not influenced by the existence of a friction. This is significant dif- 
ferent from another familiar example, the quantum Hall system, where it 
is generally understood that the energy gap in the quasiparticle excitation 
spectrum stabilizes the topological structure, meanwhile makes the dissipa- 
tion to vanish. When the quantized Hall conductance occurs it associates 
with the vanishing of longitudinal voltage. However, the total transverse 
force on a moving quantized vortex, which has a topological significance, 
remains unchanged with finite friction. 

There has been a great amount of effort in this subject, because this is 
the only type of topological singularities which can be extensively studied 
experimentally, in both classical and quantum regimes. Some of the previous 
work on vortex motion in fermionic superfluid has been discussed critically in 
[0, H. In bosonic superfluids, the situation of early work is more diverse HI. 
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So far, the experimental papers on superfluids and superconductors claim 



support to different theoretical results] 14, f5l [lq, 17]. We need to keep in 



mind that in general the real experimental situations are far more compli- 
cated than ones usually considered in theoretical models. For example, in 
the analysis of Hall measurements in superconductors, the vortex interaction 
and the pinnings are often ignored. If we use the same assumption for longi- 
tudinal resistivity, we would conclude that we should have flux flow for any 
current and at any temperature, which contradicts the longitudinal resistivity 
measurement s0. Complications also exist in mutual friction measurements. 
For example, the additional friction caused by the relative motion between 
the ends of the vortex lines and the diaphragm has not been considered in 
a recent experiment ||1 5|] . In addtion, because these vortex lines are shorter 



than the viscous penetration depth of the normal fluid, the boundary effect 
can also be important. Conclusions on the transverse force drawed from such 
experimental data should be taken with a due caution. 

Before we conclude, we discuss the thermodynamic interpretation of our 
results. We have shown that indeed a finite friction does not change the 
transverse force. To illustrate this, let us imagine a torus-shaped tank filled 
with superfluid. The normal fluid is assumed to be at rest with the tank 
initially. After creating a vortex-antivortex pair, we apply an external force 
on the vortex to let it wind around one of the circumference of the torus 
once with uniform velocity in time ttotai —* oo before their annihilation. We 
consider the work done by the external force to the system, which includes the 
superfluid, normal fluid and the tank(subtrate), and the energy and entropy 
gain of the system. The vortex is used as an intermediate step to find the 
forces. The external force on the vortex in the longitudinal direction balances 
the force felt by the vortex due to its coupling to the normal fluid. The 
external force on the vortex in the transverse direction balances the total 
force on the superfluid, which is the summation of the Mangus force and a 
possible one from the normal fluid. 

The superfluid velocity has a finite increase due to the phase change 
caused by the vortex motion. Thus the kinetic energy of the superfluid also 
increases by a finite amount. This kinetic energy increase is exactly equal to 
the work done by the external force on the superfluid, if it is the Magnus force. 
The frictions do not dissipate energy, because the process is quasi-static in 
which the normal fluid velocity is negligible. Under such a condition, the 
external force on the superfluid in the transverse direction needs to be at 
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least the magnitude of the Magnus force in order to provide the work for the 
increase of kinetic energy of the superfluid. The normal fluid and substrate 
cannot tranfer energy to the superfluid, which carries no entropy, by lowering 
their internal energies, because of the second law of thermodynamics [|18 . 



Our argument has not told us why the total transverse force is not larger 
than the Magnus force. If there is no entropy change in the normal fluid, 
then the transverse force will not be larger than the Magnus force because 
the energy increase in the system, carried by the superfluid, not normal fluid, 
must equal to the work done by the external forces. This is possible, because 
the superfluid carries zero entropy, and because our process is quasi-static, 
in which there is no entropy generation within the system. There is also no 
external heat bath to provide a heat flow. The entropy stays constant within 
the system. 

In summary, both the transverse force and friction on a moving quan- 
tized vortex has been obtained simultaneously within the framework of Ref. 
[j|]. The transverse force confirms the Berry phase results obtained with a 
trial many-body wavefunctionfQ. New results on the friction are given for a 
fermionic superfluid, which can be tested experimentally. 
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